Efficient Quantum State Estimation with Over-complete Tomography 
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It is widely accepted that the selection of measurement bases can affect the efficiency of quantum 
state estimation methods, precision of estimating an unknown state can be improved significantly 
by simply introduce a set of symmetrical measurement bases. Here we compare the efficiencies of 
estimations with different numbers of measurement bases by numerical simulation and experiment 
in optical system. The advantages of using a complete set of symmetrical measurement bases are 
illustrated more clearly. 
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I. Introduction 

Quantum states are the most important and funda- 
mental elements in the quantum world. Because of the 
superposition property, quantum states are highly com- 
plex and difficult to identify. Quantum state tomography 
(QST) [l|, [H, a universal method to reconstruct quan- 
tum states by making a series of measurements on an 
ensemble, plays an important role in quantum informa- 
tion science 0]. However, because of the asymmetri- 
cal distribution of the measurement bases of the early 
QST methods all random errors (statistical or tech- 
nological) [1] accumulate during the measurements so 
that the quantum state estimation is always inefficient 
[j| . Then a new symmetrical tomography scheme, called 
over-complete QST [(| 0] , had been established. Its ef- 
ficiency has been proven and it has been widely used in 
many experiments This method may find versa- 

tile usage in many fields, such as quantum computation 
and simulation [lj| Il 6| . quan tum communication 0, QjJ , 
quantum metrology [18|, [1911 . condensed matter physics 
20f . quantum chemistry [2ll 122^ , and quantum biology 

mm- 

In this letter, first, we analyze and numerically simu- 
late the quantum state estimation procedure and com- 
pare the efficiencies of over-complete QST with that of 
normal QST, then we experimentally demonstrate its ef- 
ficiency in optical system. 
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FIG. 1: (color online) Errors of over-complete QST and nor- 
mal QST. (a). The measurement bases of over-complete QST 
(yellow and blue) are distributed symmetrically on the sur- 
face of the Bloch Sphere, while those of normal QST (blue) 
only cover a quarter of the surface. In (b) and (c), 1000 dif- 
ferent states are chosen randomly from the entire state space. 
For each state, we simulate the tomography procedure with 
4™ (blue) and 6" (red) measurement bases for one qubit (b) 
and two qubits (c), where n is the number of qubits, with 
the same total resources N. We vary N approximately from 
10 2 to 10 6 , and calculate the error of estimation by Euclidean 
distance for each chosen value of N. Each point represents 
the mean error of estimation for the 1000 states at each dif- 
ferent N. One thousand samples are sufficient to draw the 
curve because the mean value and error bar do not change 
much between data volumes of 500 and 1000. It is obvious 
that over-complete QST is more precise. 



II. Numerical Simulation 



Suppose that there are N copies of a quantum state, 
which is totally unknown before the measurements are 
made for the reconstruction task. An n-qubit quantum 
state can be described by a 2" dimensional density ma- 
trix [25j . Because the density matrix has the properties 
of normalization, Hermiticity, and positivity, it contains 
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4™ — 1 variables, so we need at least 4™ measurements 
to determine the density matrix of the state. The stan- 
dard QST method to reconstruct the density matrix is 
to equally divide the N copies into exactly 4™ measure- 
ments, record the counts n v for every measurement, and 
then apply maximum likelihood estimation (MLE) [l[ 
to find the physical quantum state that is most likely 
to generate these measured data. To be specific, these 
measured data are assumed to be Gaussian distributed 
around the true values, so the probability of observing a 



2 



particular set of counts is given by 



P(ni, ...n m ) 



N, 



exp[-}^ — ], (1) 

lorm n aiv v 



where N norm is the normalization constant, m is the 
number of measurements, and n v — — (VvIpIVv) is 
the expectation value of \ip„) measuring the state p. 
Therefore, the physical density matrix with maximum 
P(ni, ...n m ) is most likely to be the true state. This is a 
universally applicable method to estimate any previous 
unknown quantum states. 

A serious problem of QST is that the selection of mea- 
surement bases is always asymmetric in the state space; 
in the Bloch Sphere (26| of a single qubit, for example, 
the four bases in early method cover only a quarter of the 
surface (Figure 1(a)). However, using the least possible 
number of measurement bases does not necessarily cost 
the minimum total resources. The real difficulties in ex- 
periments, such as the estimation of a quantum state of 
eight ions [27| . are the limited resources, not the measure- 
ment times. According to previous study on the choice 
of measurement sets [28j |. inscribed Platonic solid of the 
Bloch sphere measurement sets result in the minimum 
error of estimation. This indicates that symmetrically 
arranged measurements will improve the precision signif- 
icantly, in this letter, we provide insightful research on 
why the selection of measurement sets affects the pre- 
cision and experimentally demonstrate this point. By 
numerical simulation of the tomography procedure on 
single-qubit and two-qubit states, it becomes clear that 
when these N copies are distributed equally into 6 ra sym- 
metrical measurement bases, i.e., each measurement base 
receives N/6 n copies, the precision of the state estima- 
tion will be improved significantly. The chosen bases are 
of the form 



(2) 



where |W> G {|0>, |1), |+), |-), \L), \R)}, with |±) = 
^(|0> ± |1», \L) = i(|0) + \R) = i(|0) - 
which lie at the six polar symmetrical points of the Bloch 
Sphere (FIG. 1(a)). MLE is used to find the most likely 
density matrix, as in QST. The results are shown in FIG. 
1 (b) and (c), wherein the error of estimation is repre- 
sented by Euclidean distance y/Tr[(p — p) 2 ] [3j. As can 
be seen from the figure, over-complete QST needs fewer 
resources than normal QST to reach the same precision. 

In our work, the task of finding the minimum of like- 
lihood function is efficiently executed by the simulated 
annealing algorithm (29|, which proceeds as follows, (i). 
Start by estimating the target matrix, p, as an arbitrary 
physical density matrix, such as the J/4, with an initial 
temperature To. (ii). Let the current temperature T 
decrease as the procedure continues, then transform the 
matrix to another one (for instance, change an element or 
partially transpose the matrix), and calculate the value 
of the likelihood function. The probability of accepting 



the new p depends on the ratio of the new and the old 
likelihood functions, and on the current temperature T; 
more precisely, it is 



min{exp[ 



L(pn+i) - L(p n ) 
kT n 



1}, 



(3) 



where k is a constant number, (hi). Repeat step (ii) until 
the value of the likelihood function becomes stable, (iv). 
Report p as the estimated density matrix. Analogous to 
the annealing process in solid state matter, the transi- 
tion probability between two levels is determined by the 
temperature and the energy difference of the two levels. 
With T decreasing, the system will finally arrive at the 
ground state. In a similar way, the simulated annealing 
algorithm can find the minimum quickly and precisely. 

First, we analyze the situation of single qubit. Finding 
the density matrix most likely to match the experimen- 
tal data is very like the procedure of fitting a line with 
the least squares method [3(|. Obviously, at least two 
points are needed to determine a line. To determine a 
line more accurately, we may measure each of the two 
points with more resources to make the measurements 
more reliable, or we may use a series of measurements 
and fit the line with more points. The latter approach is 
better, and it is therefore always used. The reason is that 
even though any single point is not so accurate, their sta- 
tistical errors, which are random, will cancel each other. 
Moreover, if some types of random errors are related to 
the selection of the measurement bases, they will accu- 
mulate in the first method but cancel each other in the 
second. Especially when the measurement resources are 
distributed symmetrically, these errors will be expected 
to reduce to the minimum. In a similar way, when esti- 
mating an unknown quantum state, the statistical error 
is closely related to the count n v , which is determined 
by the measurement base. Technological errors, such as 
the uncertainty on the angle of the wave plate, will also 
accumulate if there are only 4 measurement bases. Sin- 
gle qubit states can be described by points on the Bloch 
Sphere, and measurement bases can only lie on the sur- 
face. Therefore, the precision of estimation will be im- 
proved significantly if we use the six polar points, which 
are distributed symmetrically on the sphere, as measure- 
ment bases. In FIG. 1(a), these measurement bases can 
be regarded as 6 directions from which to observe the 
unknown states. Any particular set of 4 of them can- 
not reach a fine estimation because they cannot see the 
state from all directions (they cover only a quarter of the 
sphere), even though they are enough to determine the 
state in theory. For more than one qubit, the space of 
states grows exponentially, and we need to use the com- 
bination of these bases described by Eq. ^ to make up 
a complete set of symmetrical measurements. 

In addition, to show more advantages of the over- 
complete QST method, we numerically simulate the es- 
timation of concurrence [311 , an important parameter in 



two-qubit entanglement verification (32|. It is defined as 



C(p) — max{0, Ai — A2 — A3 — A4}, 



(4) 



FIG. 2: (color online) Error of concurrence estimation, (a). 
Bell-diagonal state with b=0.8, the true concurrence is 0.6. 
(b). Werner state with q=0.25, the true concurrence is -0.125 
(here we do not set its minimum value to zero). N is the 
total quantity of resources used for estimation. These two 
examples illustrate that over-complete QST is always closer 
to the true value than QST. We generate the curve with 300 
samples for each point. 

where Xi are the eigenvalues, in decreasing order, of the 
Hermitian matrix R = ■Jyfppy/p-, and p is the density 
matrix. Here we do not set its minimum value to zero 
so that we may study the estimation more clearly. Two 
common states, the Bell-diagonal state [33| 

P = b\4,-){4,-\ + {i-b)\r){r\ (5) 

where \^~) = ^(|00) - |11», \4>~) = ^(|01) - |10», 
and the Werner state [34j 

p = q\ip-)(ip-\ + (l-q)I/4 (6) 

are used as examples. The results are shown in FIG. 
2. It is easy to see from the figure that the concurrence 
estimated by 36 measurements is closer to the true value 
than that estimated by 16 measurements. When data are 
insufficient, the deviation seems tremendous, the effect of 
which is caused by MLE .35]. Nonetheless, over-complete 
QST has reduced the deviation greatly. Consequently, 
the result from over-complete QST is more reliable than 
that obtained from normal QST, and over-complete QST 
is very useful to verify entanglement when resources are 
limited. 



III. Experimental Demonstration 

Now it is clear that estimation by 6™ measurements is 
better than that by 4™. Next, we experimentally demon- 
strate that a complete symmetrical set of measurement 
bases is the best choice for quantum state estimation. 
Taking two qubits as an example, here we conduct two 
all optical experiments on a Bell-diagonal state and a 
Werner state to find the optimal number of measure- 
ments for quantum state tomography. The experimen- 
tal setup is described in FIG. 3, where the polarization 
\H) and \V) are used to represent |0) and |1). Approx- 
imately 2.5 x 10 5 photon pairs are used in total, which 
are divided equally into m = 16 to m — 40 measurement 
bases for tomography (see Table 1 for the selection order 



FIG. 3: (color online) Experimental setup. Here, the abbre- 
viations of the components stand for the following: HWP - 
half wave plate, QWP - quarter wave plate, IF - interference 
filter, QP - long quartz plate, introducing complete dephas- 
ing between \H) and |V), BS - beam splitter, M - mirror, 
PBS - polarizing beamsplitter, SPD - single photon detector, 
and FC - fiber coupler. An ultraviolet doubled femtosecond 
pulsed laser (about 100 mW, 390 nm, 76 MHz) is used to 
pump two BBO crystals for type I down conversion to gen- 
erate entangled photon pairs (780 nm, -^={\HH) — IV'V))), 
and the 2nd path is split into two branches, one of them is 
controlled by optical devices to prepare either a Bell-diagonal 
state with b = 0.8 or a Werner state with q = 0.5. More 
precisely, we use an HWP to exchange \H) and |V) on the 
branch to prepare the Bell-diagonal state, while we use an 
HWP (transforming \H) and |V) to |+) and |— ), respectively) 
and QPs on the branch to prepare the Werner state. An at- 
tenuator is used to tune the ration q and b. The final state 
is detected by the automatic tomography system, which is 
controlled by a Lab VIEW program and works precisely and 
efficiently. 
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FIG. 4: (color online) Error of different numbers of measure- 
ment bases for N = 2.5 x 10 5 . (a). Werner state with q = 0.5, 
(b). Bell-diagonal state with b = 0.8. The a;-axis shows the 
number of measurement bases; the y-axis shows the Euclidean 
distance. The state estimated from 10 8 data is regarded as the 
true one. We statistically calculate the mean value and stan- 
dard deviation of the Euclidean distance from 60 individual 
experiments for each point, and we find that the mean value 
and standard deviation have already converged to certain val- 
ues. For both states, the 36-base set is the best scheme. 

of measurement bases). The results are shown in FIG. 
4. They indicate that, among the tested sets of measure- 
ment bases, the set containing 36 bases (an over-complete 
symmetrical set), is the best; the others could not achieve 
such precision because they are not distributed symmet- 
rically. This illustrates once more the importance of the 
choice of measurement bases. 

In this experiment, the quantity of total resources is re- 
garded as the sum of used resources for all measurement 
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bases. In fact, in many experiments [HI, two orthogonal 
bases, such as |0) and |1), were measured simultaneously 
because the PBS only separates them without degrading 
either one. In this way, the total amount of resources 
used by over-complete QST can be saved exponentially, 
and only 3™ measurements are required. Normal QST, on 
the other hand, with its incomplete measurement bases, 
will waste part of these resources (it also needs at least 3™ 
measurements). One caution is necessary here, however, 
saving resources in this way requires calibrating single 
photon detectors. Although a recent research [37| in- 
troduces a novel set of mutually unbiased bases for two 
qubits tomography, however, some of them are entangled 
states and require more complex setting to detect. An 
over-complete set of measurement bases in our experi- 
ment is easy to be realized in experiments and can be 
extended to more than two qubits system. 

IV. Conclusion 

In conclusion, we have made it clear that the precision 
of quantum state estimation is strongly affected by the 
selection of measurement bases. It has been proven that, 
when estimates are performed with the same quantity of 



resources and the same estimation method (MLE), over- 
complete QST makes a great difference in improving the 
precision. Recently, several improved state estimation 
methods have been proposed, such as Bayesian mean esti- 
mation (BME) 38] , hedged maximum likelihood estima- 
tion (HMLE) [231, compressed sensing method [40l |. and 
Marcus Cramer et al.'s method [4l[ for finitely correlated 
states (FCS) [11 or matrix product states (MPS) [Hj]. 
However, these methods can be only used in some spe- 
cial situations, that is, when the states are belong to some 
certain categories or prior assumption of the states are 
available. They have not revealed the intrinsic weakness 
of QST, the asymmetry distribution of the measurement 
bases. Remarkably, over-complete QST can be combined 
with these estimation methods to improve their efficiency 
significantly. We can also try other symmetrical measure- 
ment bases, such as the eight vertexes of the inscribed 
cube of the Bloch sphere. This choice of bases is similar 
to the bases of over-complete QST, and the result should 
also be very satisfactory. 
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TABLE I: CQST measurement bases selected for both experiments and simulations. Here, \L) = ^=(\H) + i\V)), \R) = 
j-{\H) - i\V)), |+) = j-(\H) + \V)), and |-> = ^(\H) - \V)). 



v measurement base v measurement base 
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\HH) 
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\HV) 
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\VH) 
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\ + R) 


10 
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11 
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13 
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\L+) 


23 


\LR) 


24 


\-R) 


25 


\HR) 


26 
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27 
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28 


\ + L) 


29 


\-H) 
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33 
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34 
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35 
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36 
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